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Abstract 

The integration of large shares of electricity produced by non-dispatchable Renewable Energy 
Sources (RES) leads to an increasingly volatile energy generation side, with temporary local 
overproduction. The application of energy storage units has the potential to use this excess 
electricity from RES efficiently and to prevent curtailment. The objective of this work is to 
calculate cost-optimal charging strategies for energy storage units used as buffers. For this 
purpose, a new mathematical optimization method is presented that is applicable to general 
storage-related problems. Due to a tremendous gain in efficiency of this method compared 
with standard solvers and proven optimality, calculations of complex problems as well as 
a high-resolution sensitivity analysis of multiple system combinations are feasible within a 
very short time. As an example technology, Power-to-Heat converters used in combination 
with thermal storage units are investigated in detail and optimal system conhgurations, 
including storage units with and without energy losses, are calculated and evaluated. The 
benehts of a problem-specihc approach are demonstrated by the mathematical simplicity of 
our approach as well as the general applicability of the proposed method. 
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1. Introduction 

Electricity generation from solar and wind resources is volatile and has limited adapt¬ 
ability to changes in electricity demand [l|. For a stable operation of the power system there 
must be a consistent balance between supply and demand. One option for achieving this 
balance is the use of energy storage units which can contribute substantially to the emansion 
of Renewable Energy Sources (RES) and their integration into existing grids 0 , 38 . 

Another possible way in which large shares of renewable energy can be integrated and 
electricity from RES can be used efficiently, especially in times of overproduction, is to 
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replace fossil fuels as the main energy source for heat by converting electricity into heat 
(Power-to-Heat, PtH) - resulting in a lower overall primary energy consumption and lower 
CO 2 emissions js], [^ . If combined with thermal storage units, PtH storage systems are a 
highly flexible option for uncoupling conversion and utilization 0,0. These heat storage 
units can be charged at considerably low losses during periods with high shares of excess 
electricity from RES and provide heat in times of low feed-in 0 . A method for storing heat 
energy used in Germany was implemented decades ago by night storage heating systems in 
private households. They operated with low overall efficiencies and subsidized electricity 
prices at hxed schedules |9[. Today, private households are still suitable for the deployment 
of PtH storage systems due to the high share of primary energy demand utilized for heating 
and hot water supply, the overall changes in the energy system and the characteristics of 


today’s heat storage units [lO 


If we assume that the electricity price represents the fluctuating availability of energy, 
these heat storage systems need to be operated in a cost-optimal manner. The optimality of 
a charging strategy for a PtH storage unit is determined by the overall electricity acquisition 
costs. The task of identifying optimal strategies is closely related to the held of mathematical 
optimization and can be described as a minimization problem. Standard solvers can be 
applied to calculate solutions, albeit with tremendously high overheads in computational 
time. Hence, the investigation of complex problems may not be possible within a reasonable 
time frame. 

The specihc structure of the minimization problem, based on a mathematical model 
of the storage units, allows the development of a new optimization method which to our 
knowledge has not yet been implemented. This work introduces an innovative optimization 
algorithm and presents proof of optimality. The solution of storage-related optimization 
problems can now be calculated in a fraction of the time required by standard algorithms. 
As a particular example, we focus on PtH storage units installed in private households. For 
these systems, cost-optimal operating strategies including storage units with thermal energy 
losses are described and used for an iterative determination of optimal system designs. 

Beyond this particular held of application, we illustrate the potential benehts of a 
problem-specihc optimization approach in terms of mathematical simplicity and compu¬ 
tational gain as compared with standard solvers. 


2. Formulation of the optimization problem 

The construction of mathematical energy storage models and the formulation of the cor¬ 
responding optimization problems, with and without constant as well as non-constant energy 
losses, are described without units in the following paragraphs. We investigate an electrical 
charging system with an attached storage unit as a buffer and consider the parameters of 
electric power consumption and storage capacity. The system needs to be connected to a 
virtual electricity grid which provides an altering price signal rehecting the availability of 
electricity. A further precondition is that the energy demand is covered exclusively by the 
system. 
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2.1. Formulation without energy losses 

Now, we formulate the mathematical model, without energy losses of the storage unit 
over time or losses during the conversion process. We discretize the considered period of 
time into n intervals of equal size. The electricity prices are denoted by Ci, C 2 ,..., c„ and 
the energy demands are represented by di, ^ 2 ,..., d„. Furthermore, we dehne the amount 
of energy used to charge the storage unit for each time-interval by xi, 0 : 2 ,..., Xn and assume 
that the prices and the demands are known. For technical reasons, the values of x are 
constrained. Each charge value cannot be negative (no discharge to the electricity grid or 
re-electri£cation) and has an upper bound C > 0, which implies: 

0<Xj<C, Vf = l,...,n. 

To cover at least the energy demand for each time-interval, additional constraints are: 

i i 

^dj Vi = l,...,n. 

i=i i=i 

The storage level for each time-interval is dehned as the difference between the quantity 
of charges and the quantity of demands up to this time-interval. Due to design-related 
restrictions, the storage level is bounded above by a maximum storage capacity value S' > 0. 
Hence, we have: 

i 

{xj — dj) < S , Vi = 1,..., n . 
i=i 

In order to minimize the total costs to cover at least the energy demand over a period of 
time of size n, we have to solve the following optimization problem: 

n 

y^^cj-xj 

i=l 

0 < Xi < C, y i = 1,... ,n (LP) 

i i i 

Xj ^ S ^ ^ dj y % — l,...,?r. 

j=i j=i j=i 

This problem is a special instance of a so-called linear program and can be solved using 
standard algorithms for general linear programs, such as the simplex method or interior 
point methods 0- 

The problem fILPp is also a special case of the problem (|P|) which is described in 
Contrary to the methods mentioned above, we utilize the special structure of (jP]) 
also of the related problem fILPp to develop a new algorithm. The basic idea of the new 
algorithm is to charge the storage unit during periods when the acquisition prices are low 
in order to avoid further purchases at times when the prices are higher. In addition to the 
price levels, the algorithm also takes into account the demand, the storage level and the 
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Appendix A 
and hence 


min 
subject to 
and 






maximum charge power for each time-interval. Therefore, the storage units are charged as 
much as possible at times of negative acquisition costs and as much as required, if the price 
is non-negative. 

The new algorithm to solve the problem ([P]) is discussed in detail in Appendix A and 
also the pseudo-code is presented (see p. [20]). Below, we present the pseudo-code of the 
new algorithm exemplary htted to the problem (ILPp . if we set Oj := di + ... + di and 
bi := S' -|- (di di) for all i = 1,..., n. If we further dehne the permutation a as 

described in Appendix A corresponding to the increasing prices by Co-(i) < ... < Co-(„), the 
pseudo-code of the new algorithm htted to the problem fILPp is given by: 


Inputs: 

ai,bi, Ci and the permutation a 

Output: 

A solution X of the problem fILPD 
for fc = 1 to n do 
Ml ■(— max {0, ai} 

i<(j(k) 

M 2 ^ max { 0 , tti} 

2>cr(fc) 

m ^ min {hi\ 

i>a{k) 

if c^(fc) > 0 then 

Xa{k) t— min I max{0, M 2 — Mi}, min{C', m — Mi}} 
else 

Xa{k) ^ minlC, m — Mi} 

end if 

for i = (j(fc), ... ,n do 
Cti ^ Cli — 
bi ^ bi 

end for 
end for 


The optimality of the new algorithm on page [12] is proven (in Appendix A ) and is one 
key element of this work. Furthermore, the source code of an implementation in Python 
is included in Appendix C This method solves problems even for large n in a fraction of 
time required by standard solvers, because at most + 3n boating point operations and 
3/2n^ -|- 5/2n comparisons to compute a solution 


see: 


Appendix Ap are required. 


To demonstrate the efficiency of the new algorithm, we compare the runtimes between 
the algorithm and the common solver linprog in the following. For this comparison, a 
straightforward implementation of the new algorithm in Python (cf. Appendix C ) and the 
linprog implementation as available in MATLAB 2015b were used. The calculations were 
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n 

100 

500 

1000 

2500 

5000 

Runtime algorithm [s] 

0.006 

0.03 

0.07 

0.18 

0.4 

Runtime linprog [s] 

1.008 

2.41 

18.16 

232.31 

1308.9 


Table 1: Calculation times of the new algorithm compared with linprog (MATLAB) of the 
same problem fILPjl . with respect to the problem size n. In addition to the general efficiency 
of the new algorithm, its benehcial scaling behavior (runtime increases by a factor 65 if the 
problem size is increased from 100 to 5000, whereas the scaling factor for linprog is about 
1300) becomes evident. 


performed on a desktop computeili], based on input data for the problem ([P]) as described 
in section 13.11 The results are presented in table [T] and emphasize the efficiency of the new 
algorithm and its benehcial scaling with the size of the problem. This efficiency among 
others allows to perform detailed sensitivity analyses and solve highly complex optimization 
problems, where problem ([P]) has to be solved thousands of times as a sub-problem (cf. 
section ED. For instance, the calculations performed in section 13.21 would have taken more 
than 60 days by using the solver linprog compared with about 25 minutes for the newly 
proposed algorithm. 

To include energy losses of the storage unit or during the conversion process, which can 
be constant or not-constant, the constraints of flLPp have to be transformed into the shape 
of the general problem ([PD, before the new algorithm can be applied. These modihcations 
are presented below for the two different types which can also be combined. 


2.2. Constant energy losses 

To include constant energy losses h,... ,ln for each time-interval, the constraints in fILPp 
have to be modihed. These losses can simply be added to the energy demand for each 
time-interval, which leads to the following problem, with the same structure as ( 0 : 

n 

min E Ci ■ Xi 

i=l 

subject to 0 < Xi < C , Wi = 1,... ,n (1) 

i i i 

and dj + li < Xj < S + dj + li^ Vi = 1,..., n. 

j=i j=i j=i 


2.3. Non-constant energy losses 


A slightly more complicated modihcation of flLPp involves the integration of non-constant 
energy losses, which is described in Appendix A.l in detail. This represents for instance 


1 Details: Intel® Core^“ i7-930 Processor (2.80 GHz), 12GB RAM, MATLAB 2015b, Python 2.7. 
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losses during the charging process or energy losses of the storage unit over time. The modihed 
problem is: 


n 


min 

y^^cj- Xj 

2=1 



subject to 

0 < Xi < C, Vi = l,.. 

.,n 


and 

i i 

Y^q'-> ■ d, <Y,q<-> 

i 

■ Xj < S + ■ dj , Vi = 1,.. 

.,n 


j=i j=i j=i 


( 2 ) 


with an energy loss factor g > 0. This problem can be transformed into an equivalent 
problem of the form ([P]), which is described in Appendix A.l in detail. 

Note that the minimization problems flLPp . ([T]) and (|2]) are special instances of the 
problem (jP]) (paged!]). Because the new algorithm (pagedOj) solves the general problem (iPj) . 
it also solves the above mentioned problems. Thus, this algorithm will be used for all 
discussions and calculations in the following investigations. 


3. Applications to PtH storage units and numerical results 


The technical realization of a PtH conversion unit can be achieved for example by an 
electric heat pump system or an electrical heating element, using the concept of an immer¬ 
sion heater. Common conversion systems offer almost 100 % efficiency for a simple heating 
element or even more for a heat pump system, with a performance factor of three to four 


12, 13 


Therefore, we do not include energy losses of the PtH system during the con¬ 
version process. The thermal storage unit, which is assumed to be attached to the PtH 
system, can be realized for ex amp le by one of these three storage groups: sensitive, la¬ 


tent, or thermochemical BQ Il5| . For simplicity, the relevant parameters are the electric 
power consumption for the converter, the different storage capacities and the storage unit 
efficiencies. 

This section presents a procedure to determine an example of a heat demand prohle 
and calculations of cost-optimal charge curves via the new algorithm for storage units with 
and without energy losses. Furthermore, we present an approach to estimate cost-optimal 
combinations of heating systems and connected storage units, with different efficiencies of 
the storage unit. 

To perform a detailed analysis of the demand for thermal energy in private households, 
more accurate values than the specihc heat demand for one year are required. A common 
option for calculating these load prohles is to use the VDI Guideline 4655 “reference load 
prohles of single and multi-family homes for the use of CHP plants”. The guideline includes 
the following influence variables to determine such prohles: type of household, number of 
residents, annual energy demand for space heating and hot water supply as well as the 
climate zone. 

The majority of households in Germ any in 2011 were single and double-family homes, 
with an average living area of 116.6 m^ 16|. The number of households with two residents 


6 













3 



0 -^^^^^^^^^^^- 

Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul 


Figure 1: Hourly heat demand (space and water heating) profile for a single-family household 
with two persons, living space of 100 m^, specific heat demand of 100 kWh/(m^-a) and annual 
heat demand of 10741.15 kWh/a, (01.07.2012 - 30.06.2013). 


ranked second with more than 33 
for the future 


171,118 


in the same year and with the highest rate of increase 
In this work, we study a single-family building with two residents. 


The specific heat demand of these households was on average around 150 kWh/(m^-a 


m 


2011, but displayed a decreasing tendency when compare to previous years [19|]. With 


respect to the EnEV 2009 and the sustained decreasing specific heat demands, a value of 
100 kWh/(m^-a) is considered (l^. 20 |. The guideline recommends an annual heat demand 


for hot water supply of 500 kWh/(person-a), which depends considerably on the number of 
residents of the observed household and their habits. However, it can be assumed that this 
value underestimates today’s consumption and the adoption of a signihcantly higher number 
for real averages approximates this value more closely. In this work, we assume a value of 
750 kWh/(person-a). 

Furthermore, we set the observed period of one year from July 2012 to July 2013 and used 
the records of the German Meteorological Service for the weather station Bremen (airport 
BRE) as the data basis for the climate conditions to assign each day with a certain type of 
weighting factor, contrary to what is envisaged in the VDI 4655. The annual heat demand 
for the chosen household is 10741.15 kWh/a; this is about 6 % less than the theoretical 
guideline value. For a higher resolution, the total demand values for each day are further 
divided equally to observe the hourly demands depending on the daily energy requirements. 
From the preliminary considerations we calculated an example load profile for the following 
household: single-family household (short: SFH) with 2 persons, living space of 100 m^, 
specific heat demand of 100 kWh/(m^-a) and an annual heat demand of 10741.15 kWh/a, 
which is illustrated in figure [U 

As a data basis for the electricity acquisition costs representing the volatility of the RES 
feed-in, we used the Physical Electricity Index (PHELIX). This reference price signal for the 
observed period published on the website http://www.epexspot.com for auction results of the 
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Figure 2: Hourly price development of the electricity price (PHELIX) at the spot market 
EPEX SPOT, (01.07.2012 - 30.06.2013). 


PHELIX and is illustrated in figure [2] with an average price value of about 0.04 EUR/kWh. 
3.1. Cost-optimal charging strategies 

We consider the SFH for two different PtH storage systems in the following. To evaluate 
and compare the examples more efficiently and even better, we set the maximal storage 
capacities for the thermal storage units as relative numbers corresponding to the mean daily 
average heat demand. The different observed systems are first a common PtH storage system 
consisting of a half-day storage unit (~ 14.71 kWh) and a PtH conversion unit with a maxi¬ 
mal charge power of 9 kW and second a system defined by a 14-day storage (~ 411.99 kWh) 
and a conversion unit with a maximal charge power of 25 kW. This system neglects the 
switching limits of today’s power connections and the today’s common dimension of heat 
storage units in private households, but is chosen to demonstrate the cost saving potentials 
and also the overall benefit of a cost-optimal operation in combination with systems be¬ 
yond today’s standards. As mentioned above, the efficiencies of both conversion units for 
the moment are set to 100 %. These observed systems are denoted as SFH 1/2 and SFH 14 
respectively. As a reference value for the annual electricity acquisition costs and in order 
to compare the different cost savings, the costs with no storage unit would be 443.13 EUR 
to cover the annual heat demand; this represents the conversion of electricity into heat at 
times when the energy is required. 

The calculated curves in figure [3] for SFH 1/2 and in figure 0] for SFH 14 result from the 
new algorithm and describe the cost-optimal strategy for operating the modeled PtH storage 
unit without energy losses. For SFH 1 / 2 , the overall annual electricity acquisition costs are 
310.99 EUR, which is 29.82 % less than the annual costs resulting from a system with no 
storage unit. The overall costs for SFH 14 are 71.55 EUR, which is more than 83 % below 
the reference value. The increased heat demand between November and April (cf. figure 0]) 
is represented by a higher number of charges within this period to cover the demand and 
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by a highly fluctuating storage level. This behaviour is much more pronounced in the result 
of SFHi /2 compared with SFH14, which is due to the smaller storage capacity and also the 
lower charge power (cf. figures [3] and 0]) . As described, the storage unit is charged as much 
as possible during hours with negative prices. This leads to a storage level greater or equal 
to zero at the end of the observed period. The storage level of SFH 1/2 is zero at the end of 
the period (cf. figure [3]), compared to a storage level of 146.32 kWh for SFH 14 (cf. figure 
i). As a consequence, the large PtH storage system utilizes the negative electricity prices 
considerably better than the common smaller system. 

The difference between both storage capacities is also reflected in the decreasing number 
of non-zero charges from 2872 for SFH 1/2 down to 514 for SFH 14 as the capacity and charge 
power increases. For SFH 1 / 2 , the total number of charges with maximum charge power is 
575 and more in comparison to 407 maximum charges for SFH 14 (cf. figure |5]). The number 
of charges, which are non-zero and less than the maximum charge power, turn out to be 
very different (cf. figure |5]). These are 2297 for SFH 1/2 compared with 107 (less than 5 %) 
for SFH 14 . In summary, the storage capacity and the charge power have a strong impact 
on the behaviour of the resulting cost-optimal charge curve for the chosen household. As 
these factors rises, the number of charges lower than the maximum charge power decreases 
drastically and the number of charges with the maximum power increases slightly. This 
leads to a shift in the values of the cost-optimal charge strategy towards the maximum 
charge power. This fact is presented in figure [5a] for SFH 1/2 and in figure |5b] for SFH 14 . 
Further, increasing the system factors reduces operating hours tremendously, which leads in 
combination with a storage level higher than zero at the end to a sharply higher cost saving 
potential for large systems. 

We now add energy losses of the thermal storage unit to the problem as described in 
section |2] and further maintain no energy losses of the conversion system. With respect to 
the directive 2010/30/EU of the European Parliament and of the Council j^, we applied 


an exponential regression analysis on the data taken from j 22 |, to calculate heat losses in 


standby mode (details in Appendix B ). The result of the regression leads to an hourly factor 
of about 0.9962 (approx. 8.73 % daily energy losses) for the half-day storage unit and an 
hourly factor of about 0.9989 (approx. 2.61 % daily energy losses) for the 14-day storage 
unit. 

The differences between the cost-optimal charge curve without energy losses and the 
solution of the problem with losses for SFH 1/2 are quite small. The charges of the different 
charge curves are shifted by only a few hours in time, which means that the moment in time 
and the amount of energy are almost identical for both examples. As a result, almost no 
differences between the cost-optimal charge curves of a heat storage system with a lossless 
and a lossy small storage capacity. By contrast, the differences in the SFH 14 versions display 
a block-like structure. The behaviour is either maximum charge power for the cost-optimal 
charges with a lossless storage unit or for the lossy version over a couple of weeks (cf. 
figure [ 6 ]). The annual electricity acquisition costs for SFH 1/2 with a lossy storage unit is 
318.09 EUR and 96.78 EUR for the lossy version SFH 14 . This means an increase of about 
2.28 % for the first and 35.26 % for the second system. Nevertheless, it is 28.22 % and 
78.16 % respectively less than the reference value. This is caused by the additional amount 
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Figure 3: Cost-optimal hourly charge curve (top panel) and corresponding hourly storage 
level (bottom panel) for a storage unit with capacity 14.71 kWh, (01.07.2012 - 30.06.2013). 
The solution is 310.99 EUR for the electricity acquisition costs. (Assumptions: single-family 
household with 2 residents, living area of 100 m^, specific heat demand of 100 kWh/(m^-a), 
storage capacity of 14.71 kWh (1/2-day of average heat demand) and charge power of 9 kW, 
no thermal energy losses) 


of energy needed to cover the heat demand, which is 200.78 kWh for SFH 1/2 and 994.19 kWh 
for SFH 14 (cf. hgurelHD- 

Below we address the question of how the efficiency factor affects the solution. For the 
two example systems with a charge power of 9 kW and 25 kW respectively and two different 
storage capacities of half a day and 14 days as well as the household specifications as above, 
the results of a variable hourly efficiency of the thermal storage unit for both examples are 
presented in figure[71 The hourly efficiencies are between 0.95 and 1.00, which represent daily 
energy losses of between 70 % and 0 %. The annual electricity acquisition costs decrease as 
the efficiency increases in both cases. This is due to the lower amount of energy required for 
an improved isolated storage unit. The curvature of the curve for the 14-day storage unit 
is growing faster as the efficiency increases, in comparison with the curve for the half-day 
storage unit. This leads to higher influences of the efficiencies for larger storage units, which 
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Figure 4: Cost-optimal hourly charge curve (top panel) and corresponding hourly storage 
level (bottom panel) for a storage unit with capacity 411.99 kWh, (01.07.2012 - 30.06.2013). 
The solution is 71.55 EUR for the electricity acquisition costs. (Assumptions: single-family 
household with 2 residents, living area of 100 m^, specific heat demand of 100 kWh/(m^-a), 
storage capacity of 411.99 kWh (14-day of average heat demand) and charge power of 25 kW, 
no thermal energy losses) 
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Figure 5: (a) Histogram of the cost-optimal charges with a maximum charge power of 9 kW 
and storage capacity of 14.71 kWh (1/2-day) and an efficiency of 100 %. (b) Histogram of 
the cost-optimal charges with a maximum charge power of 25 kW and storage capacity of 
411.99 kWh (14-day) and an efficiency of 100 %. 



Figure 6: Differences between both cost-optimal charging strategies for a storage unit with a 
capacity of 411.99 kWh (14-day) and hourly efficiencies of 100 % and 99.89 %, respectively. 
Positive values indicate higher charges at the respective times in the 100 %-case. As a general 
trend decreased efficiencies result in increased energy consumption and delayed charging of 
the storage. The additional annual amount of energy required for cost-optimal charging of 
the storage unit with hnite efficiency of 99.89 % is 994.19 kW. 
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Hourly efficiency [%] 

Figure 7: Electricity acquisition costs for PtH storage systems with a maximum charge 
power of 9 kW and a storage capacity of half a day (~ 14.71 kWh, dashed orange line) and 
also 25 kW charge power in combination with a storage capacity of 14 days (~ 411.99 kWh, 
solid blue line), both for hourly storage efficiencies from 95 % up to 100 %. The calculated 
factors of energy losses for these storage capacities are about 99.62 % (red triangle) for the 
1/2-day storage unit, which leads to overall costs of 318.09 EUR and 99.89 % (red cross) for 
the 14-day storage unit, which results in overall costs of 96.78 EUR. 


have been indicated in the previous examples. In the case of hourly efficiencies above 0.99 
(daily energy losses of about 22 %) in particular (cf. £gure[7]), the optimal annual costs for 
larger capacities decrease more quickly compared with the trajectory of small ones. 

3.2. Operating costs for different system combinations 

Now we present the application of sensitivity analysis to the optimization problem and 
point out the benehts of the new algorithm. The considerations in 13.11 were based on given 
(hxed) values C and S for the technical restrictions on the conversion unit and storage 
capacity. In addition to the acquisition costs, installation costs for a PtH storage unit will 
be observed from now on. The aim is to minimize the annual electricity acquisition costs, 
with the inclusion of costs for the purchase and installation of different PtH storage units. 
This more economic approach for estimations of the overall costs leads to the minimization 
problem: 


min /C (C) + £ (5) + M (C, S) , (3) 

C,S<oo 

with (non-linear) functions /C(C), £(*5) and M(C,S) as the initial minimization prob¬ 
lem flLPD . Here, /C (C) represents the cost-function of a PtH unit depending on a certain 
maximum charge power C and £ (iS) the cost-function for a storage unit with capacity S. 
The solution of ([2]) requires C times S calculations of the optimization problem ([P]). Due 
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storage capacity [kWh] 

Figure 8: Total costs for different system setups. The minimum value is 280.02 EUR for 
a system with C = 28 kW and 5 = 70 kWh (about two and a half days). The 5 % level- 
line (red dashed line) is at 294.02 EUR and the 10 % level-line (red solid-dashed line) at 
308.02 EUR. (Assumptions: life span of 15 years, charge power between 3 kW and 100 kW, 
storage capacity between 0 kWh (no storage unit) and 410 kWh (14-day storage unit) and 
an efficiency of 100 %) 


to the gain in efficiency of the new algorithm (see: Appendix A), solving the problem fILPp 
for even thousands of different combinations C, S is possible within short time frame (cf. 
tabled]). 

The household data set, the observed time-period and the price signal are the same as in 
[Q and the life span of the installed PtH storage system is assumed to be 15 years. For the 
time being, the efficiencies of the storage units are set to 100 %. For ease of presentation 
maintenance costs and sales taxes are not included in the analysis at this stage. 

Based on regressions detailed in [Appendix B[ the partial cost-functions K. (C) and £ (S) 
for the annual costs of charging and storage units can be approximated from data, resulting 
in the specihc optimization problem 


min (0.47 ■ C EUR/kW) + (0.95 ■ 5 EUR/kWh) + M (C, S) . 


(4) 


The maximum hourly demand of the example household is about 3 kWh, so we set 
3 < C < 100, where the upper bound may not be feasible in practice and only repre¬ 
sents a theoretical value to analyze the behaviour of the solutions. The observed storage 
capacities are between zero, which represents a system without any storage unit, and a stor¬ 
age capacity of 410 kWh, which is about the 14-day daily mean average heat demand of the 
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household. The subdivision of S is set to 10 kWh 8 hours) and the separation of C is set 
to 1 kW. 

The contours of the solutions for each combination are illustrated in figure [HI The 
solution of the minimization problem (jl]) is 280.02 EUR for the combination C = 28 kW 
and iS = 70 kWh (cf. figure [H]). It has to be mentioned that this solution does not represent 
additional installation, operational and further fix costs as discussed above. In terms of 
sensitivity of the solution to the problem (jl]), we study the area surrounding the minimum. 
The 5 % and 10% levels, which represent the area of solutions that are 5 % and 10% 
respectively higher than the optimum, are 294.02 EUR and 308.02 EUR (cf. fig. |8]). This 
comprises a wide range of the calculated points and is about 15.74 % of all calculated points 
for the 5 % level as well as more than 30 % for the 10 % level. These areas, which represent 
about the same overall costs, can be used to estimate economic combinations for PtH systems 
and connected storage units. A result of this approach is that the combination for a charge 
power of 6 kW and a storage capacity of 30 kWh results in about the same total costs for 
one year as a combination of 100 kW for the charge power and 160 kWh for the storage 
capacity. This leads to a relatively robust and insensitive behaviour of the solution, but 
depends on the functions fC (C) and £ (S) in (jl]). Further, hgure [8] shows that the inclusion 
of installation costs reduces the cost savings potentials especially for large storage systems, 
in contrast to the results in 13.11 


3.3. Operating costs for different system combinations including losses 

In 13.21 we considered the storage units with an efficiency of 100 %. We now add energy 
losses of the storage unit to the problem (jl]), which will be dependent on the maximum 
storage capacity. Therefore, the function A4 (C,S) in problem (jl]) is now represented by ([2]). 
We use the same procedure as described in Appendix B to calculate hourly energy losses for 
each storage capacity. The data sets are assumed to be the same as in l3.2l but with additional 
loss factors for the different storage capacities 5 G {0 kWh, 10 kWh, 20 kWh,..., 410 kWh}. 
The contours of the solutions are illustrated in hgure |9l 

The solution of the minimization problem with lossy storage units is 291.50 EUR for the 
combination C = 26 kW and iS = 60 kWh (cf. hgure [9]), which is just 4.1 % higher than 
the solution in 13.21 Again, we study the sensitivity and robustness of the solution. The 
5 % level has a value of about 306.08 EUR and the 10 % level is 320.65 EUR. As a result, 
the consideration of storage units with energy losses leads to a solution which is close to 
the solution of (jl]), but the area with equal solutions for diherent combinations is smaller. 
Therefore, the overall costs for combinations with energy losses increase faster away from 
the minimum than the costs without losses, but are still relatively insensitive. 


4. Conclusions 

The integration of high shares of electricity generated by RES into existing electricity 
grids constitutes a major challenge for future energy systems. Energy storage units can 
contribute signihcantly to the expansion and integration of electricity from RES. 
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Figure 9: Total costs for different systems with the same representation and assumptions 
as described already in figure [HI but now with thermal storage losses explicitly included in 
the analysis (see section for details). The resulting cost-minimum is 291.50 EUR for a 
system with C = 26 kW and 5 = 60 kWh (about two days). The 5 % level-line (red dashed 
line) is at 306.08 EUR and the 10 % level-line (red solid-dashed line) at 320.65 EUR. 


In section O we developed a mathematical model to describe conversion systems com¬ 
bined with storage units. Based on the idea of using electricity from RES efficiently, utilizing 
local overproduction and increasing the share of RES within electricity grids, we assumed 
that the electricity price represents the availability of electrical energy as well as the volatile 
feed-in from RES. Next, we investigated the cost-optimal charging of these energy storage 
units, which is linked to the minimum electricity acquisition costs. This minimization prob¬ 
lem was expressed as a mathematical optimization problem. Furthermore, a procedure for 
the inclusion of constant and non-constant energy losses was described, which covers a wide 
range of applications. 

The focus of this work was the development of a new optimization method to calculate 
solutions for the modeled storage-related optimization problems and proof of optimality. 
Due to the special structure of the storage problem, this method provides a tremendous 
gain in efficiency and is also applicable to general storage-related optimization problems. 

In section [31 we applied the approach from the previous section to the PtH technology 
combined with thermal storage units. The aim was to study the principal effects of our 
approach and the consequences on the resulting cost-optimal charging strategies. Based on 
a fluctuating electricity price as a measure for the electricity acquisition costs, we calculated 


16 




















cost-optimal charging strategies for different PtH storage systems including storage units 
with and without energy losses. We found that the increase in capacity leads to a clustering 
of the cost-optimal charge values. Most of them are either zero or maximum charge power, 
which leads to fewer fluctuations of the storage level over the observed time period. Fur¬ 
thermore, we evaluated the impact of a variable efficiency factor on the minimum electricity 
acquisition costs for two different storage capacities, resulting in a higher efficiency impact 
for large storage units. We also described an approach to estimate cost-optimal combina¬ 
tions of the conversion system and the storage unit (with and without energy losses), with 
the inclusion of installation costs. The result is that the inclusion of energy losses decreases 
the cost-optimal solution only by a small amount and retains the elliptic shape compared 
with the lossless solutions. This leads to limited influences of different efficiencies on the 
cost-optimal combination for the observed household. But the inclusion of installation costs 
signihcantly decreases the cost savings potentials especially for large PtH storage systems. 
At the end of this section, we illustrated the robustness of the overall costs with respect to 
the optimal solution, as a consequence of the applied sensitivity analysis and have shown 
that this optimal solution can be reached by systems, which are larger than today. 

In this work, we used the fluctuating electricity spot market price as a measure for the 
approximation of the electricity acquisition costs. In practice, this measure underestimates 
common acquisition costs. On the contrary, the utilization of a conversion process with 
a high performance factor, e.g. a heat pump system, could have a compensatory effect. 
Nevertheless, the investigations in this work illustrate the resultant effects of a cost-optimal 
charging strategy. 

In conclusion, the novel optimization method developed in this work is applicable to 
storage-related optimization problems, even where different efficiencies are included. Due to 
the reduction in runtime, the presented method can also be used to manage a collection of 
decentralized storage units, for instance via the utilization of small ECUs. Furthermore, this 
work reveals the benehts of a problem-specihc approach in terms of mathematical simplicity, 
applicability to general storage-related problems and decrease of computational effort for 
even high-resolution problems on a large time-scale. 
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Appendix A. New algorithm: definition and proof of optimality 


Let Oj, bi, Ci,Ui G M with m* > 0 and Oj < bi for i E I n}. We are interested in 

finding a point x = (xi,... ,x„) G M” which solves the following constrained minimization 
problem flPj) . 


min E Cj ■ Xj 
snbject to 0 < x* < 

i 

and tti < Xj < bi for alH G /. 


(P) 


We denote its minimnm valne by Pmin- A feasible point of dP]) is any x G M” which fnlfills 
the constraints, and a solution of ([P]) is a feasible point x which achieves the minimnm valne 
Pmin, i.e. ■ % = -Pmin- Note that, althongh Pmin is nniqne, there conld be different 

solntions x if some of the q coincide or are eqnal to zero. We assnme that there exists 
at least one feasible point. Hence (jP]) has at least one solntion, because the constraints 
0 < Xj < Mj ensure that the set of all feasible points is compact. For more information about 


the solvability of linear programs we refer to [ll|. There the reader can also find common 


algorithms to solve linear programs, where among the best known are the Simplex Method 
and Interior Point Methods. Here we exploit the special structure of (jP]) to develop the 
much more efficient new algorithm, which needs at most + 3n floating point operations 
and ?)/2n^ + 5/2n comparisons to compute a solution. The source code of this algorithm 
written in the programming language Python is presented in [Appendix C 


We now give a precise description of the new algorithm, which is denoted as Algorithm 
in the following. Let a be a permutation of the indices which correspond to increasing Cj 
values, i.e. Co-(i) < ... < Co-(„,), and let x be a solution of ([P]). Then Algorithm successively 
computes the values Xa-(i),..., Xa{n)- In fact we show in the following 5 steps that the value 
Xcr(fc) must necessarily be equal to the explicitly computable value in Line 6 or 8 of Algorithm, 
which proves that Algorithm indeed computes a solution of ([P]). At first we assume that all 
Cj are different, i.e. Co-(i) < ... < Co-(n), and not equal to zero. 

Step 1 ; We show that Pmin = Co-(i) • Xo-(i) + Pmim where Is the minimum value of the 
reduced linear program 


min 
subject to 


Co-Q) • 

j=2 

0 < Xi < Ui 


and a[ < Xj < b[ for alH G / \ {cr(l)} , 

jVo-(I) 


(P’) 


19 





Algorithm 

Input: Oj, bi, Ci, Ui and the permutation a 
Output: A solution x of the problem ([P]) 

1 ; for k = 1 io n do 
2 : Ml ^ max{ 0 ,aj} 

i<(T(k) 

3: M 2 ^ max{0,ai| 

i>a(k) 

4: m ^ min {hi\ 

5: if C^(fc) > 0 then 

6 : Xa{k) ^ min | max{ 0 , M 2 — Mi}, min{no.(A:), m — Mi}} 

7: else 

8: Xa{k) ^ min{M^(fc), m - Mi} 

9; end if 

10 : for i = a{k),... ,n do 

11 . di i Q/i Xfj(^k) 

12: bi^bi- 

13: end for 

14: end for 


with the modihed lower and upper bounds 


a 


/ 

i 


b 


/ 

i 


{di 

max}ci^p^_i, } 

di ^cr(l) 

(h 

< min{6^(i)_i,6^(i) -£^( 1 )} 
Xa{l) 


, i < (t(1 ) — 1 
,i = a(l) - 1 
,i > a(l) 

, i < cr(l) — 1 
,i = ct( 1 ) - 1 
> ct( 1 ) 


Indeed, the point (hi,..., ho-(i)-i, ha-(i)+i,..., Xn) G M” ^ is obviously feasible for flP^ . and 
hence we have 


Pmin ^cr(l) ' h(j(i) 4“ ^ ^ — ^o'(l) ' ^o'(l) ~k Pfrii', 

i=2 


Conversely let (xi, ..., a;o-(i)-i, a;o-(i)+i, ..., Xn) G M” ^ be a solution of flP^ . Then the point 
(xi,..., Xo-(i)-i, ho-(i), Xo-(i)+i,..., Xn) G M” is feasible for (jP]), and hence we also have 


Co-(l) ■ h(j(l) T Pmin ' ^cr(l) T ^ ^ ^cr(i) ' — Pmin ■ 

i=2 


From this we moreover infer that (hi,..., ho-(i)_i, ho-(i)+i, • • •, h„) is actually a solution 
of flP^ . As a consequence, once the value of ho-(i) is known, it remains to compute the 
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remaining values Xcr( 2 ), ■ ■ ■, Xo-(n) by solving flP^ . Since has the same structure as o. 
we can proceed analoguously to get -P^m = ^o-( 2 ) ■ £o-( 2 ) + Pmin^ with a similar problem (P”), 
and determine the value of £o-( 2 )- Iterating this process n times leads to the outer for-loop 
in Algorithm. The inner for-loop corresponds to the modification of the lower and upper 
bounds in flP^ (note that only maxima of the and minima of the bi are needed to compute 
Xa-(k))- If remains to show that fo-(i) equals the explicit value in Line 6 or 8 of Algorithm, 
because then this inductively holds for all Xa{k)- As in Line 2-4 of Algorithm we define 

Ml := max{ 0 ,aj} , M 2 := max{ 0 ,aj} , m := min {feA . 

^<(t( 1 ) 

Step 2: We show that ^ Afi and 0 < f (t(i) < min{Mo-(i),m — Mi}. The hrst 

inequality holds because for alH < (j(l) we have 


i 

xj>Y % ^ • 

i=i i=i 


And the second inequality holds because of 0 < x„(i) < Mo-(i), and because for all i > (j(l) 
we have 

i 

h>Y • 

i=i i=i 

Step 3 : We show that in case Xa-{i) < min{Mo-(i),m — Mi} we must have XljiV ^ ~ Afi. 

Assume to the contrary that % > Afi > 0. Then we find an index i < (t( 1) such that 

Xi > 0. Let ii be maximal with this property, i.e. Xi^ > 0 and Xi = 0 for all ii < i < ct( 1). 
Then we hnd some 5 > 0 such that — 5 > 0, Xj — 6 > Mi and Xo-(i) + 5 < Mo-(i)- 

Now we define x G M"" by 

{ Xji S , i ii 

£a(i)+(5 A = cr(l) • 

Xi , otherwise 


Then for all i < ii and all i > (t( 1) we have And for all ii < i < ct(1) 

we have 


i i 

j=i j=i 

Hence x is feasible for ([P]) with 



< bi 

= Xj - S > Mi> Qi 



Xj Pmin T b • (Cct(1) 


Qi) ^ Prr 


which is a contradiction to the fact that f is a solution of flPj) . 

Step 4: In case Co-(i) > 0 we prove that 

Xo-(i) = min { max{0, M 2 — Mi}, min{Mo-(fc), nr — Mi}} . 

For this we consider three cases: 
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1. In case Xa{i) = M2 — Mi it follows with the second inequality in Step 2 that 0 < 
M2 — Ml < niin{Mo-(i), m — Mi}, and hence the assertion is true. 

2. In case a;o-(i) < M2 — Mi we show that Xcj{i) = niin{Mo-(i),m — Mi}, from which the 

assertion follows. Assume to the contrary that Xa(i) < min{Mo-(i),m — Mi}. Then 
it follows from Step 3 that Let Iq > (j(l) be an index with 0*0 = 

maXj>o-(i){aj}. Then we have 


iQ iQ iQ 

M2 < ^ Xj = Ml + £<^(1) + ^ Xj < M2 + ^ Xj , 

J = 1 j = o-(l) + l j = CT(l) + l 


from which we infer that X]7=o-(i)+i % ^ H^iice we hnd an index i > cr(l) such that 
Xi > 0. Let ii be minimal with this property, i.e. hjj > 0 and f j = 0 for all (t(1) < i < 
ii. Then we hnd some <5 > 0 such that f — 5 > 0 and fo-(i) + ^ < min{Mo-(i), m — Mi}. 
Let a: G M” be dehned as in Step 3. Then for all i < (t( 1) and all i > A we have 
Xj. And for all (t(1) < i < A we have 


i 


i 


= 

i=i 




i=i 


> tti 

= Ml + x„(^i) + 5 < m <hi 


Hence x is feasible for ([P]) . This leads to the same contradiction as at the end of Step 

3. 

3. In case 0 :^( 1 ) > M 2 — Mi we show that x„{i) = 0, from which the assertion follows. 
Assume to the contrary that x„(^i) > 0. Then we hnd some 5 > 0 such that x„{i) — 5 > 0 
and Xo-(i) — 5 > M2 — Mi. Now we dehne a: G M" by 


Xi 


Xa{i)-5 ,i = a{l) 
Xi , otherwise 


Then for all i < cr(l) we have %• for all i > cr(l) we have 


i i 

j=i j=i 

Hence x is feasible for ([P]) with 



< h 

> Ml + ai(j(i) — S > M2 > cii 



P - S ■ c 

^ m/i.n. ^ 


o-(l) 


i&i 


< Pn 


which is a contradiction to the fact that a; is a solution of flPj) . 

Step 5: In case Co-(i) < 0 we prove that fo-(i) = niin{MCT(fc), m —Mi}. Assume to the contrary 
that Xa{i) < min{Mo-(i),m —Ml}. Then it follows from Step 3 that If there 
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is some index ii > (t( 1) such that Xi^ > 0, then we can argue as in the second case of Step 4 
to arrive at a contradiction. Hence we have h* = 0 for all i > cr(l). Let 5 > 0 be such that 
Xa(i) + 5 < min{Mo-(i),m — Mi} and dehne a; G M"" by 


Xi := 


£< 7 ( 1 ) +5 ,f = cr(l) 
Xi , otherwise 


Then for all i < cr(l) we have * — *^(1) have 


i=i i=i 

Hence x is feasible for ([P]) with 


> tti 

= Ml + Xa{i) + 5 < m <bi 


/ ^ Cj • Xj Pmin T ^ • C(j(l) ^ Pmin i 
j&I 

which is a contradiction to the fact that x is a solution of flPjl . 

The above proof reveals that, when the q are different and not equal to zero, the solution 
X of ([P]) is in fact unique and does not depend on the actual values Cj but only on their 
order via a and whether they are positive or negative. Note that Algorithm only uses this 
information, too. It remains to show that Algorithm also computes a solution of ([P]) if some 
of the Ci coincide or are equal to zero. Fix a permutation a of the indices which corresponds 
to increasing Cj values, i.e. Co-(i) < ... < Co-(„), and dehne := Co-(j) + j ■ 5 for 5 > 0 and 
all j G /. Then we still have the same, but now strict, order < ... < and for 

all 5 > 0 small enough we also have ^ 0 for all j E I and positive (negative) values 
remain positive (negative). Hence for all 5 > 0 small enough all linear programs (P*^) with 
Ci replaced by cf have the same unique solution x, independent of 6. We show that x is also 
a solution of ([P]). Indeed, since o and (P^) have the same constraints, for any solution x 
of dE]) we have 


Pmin = ^ - Xj c] -Xj-S-^J-Xj. 

jei jei j£i jei jei 

By letting 5 —)■ 0 we see that the right hand side converges to ' Xj = Pmin- Hence we 

also have ' Xj = Pmin, be. x is a solution of ([P]). Note that in this way a zero value 

Cj = 0 is treated in the same way as a positive value, which is also rehected by the choice “>” 
in Line 5 of Algorithm. Based on what we have said above there is no need to change the 
input values from Cj to cf, since the actual values do not matter. Furthermore, we could as 
well use “>” in Line 5 of Algorithm to treat zero values in the same way as negative values, 
and a similar proof as above with := Co-(j) — | conhrms that this modihed Algorithm 
also computes a solution of flPl) . 
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Appendix A.l. Modifications 

A slight modification of ([P]) with constraints like 

i 

pj ■ Xj <bi, W i e I, 
i=i 

and positive numbers e,p > 0 can be transformed in an equivalent linear program of the 
form (P), by setting 

~ ~ Cj ^ Oi ~ bi ^ 

Xj := Pj ■ Xj , Cj := — , Uj := — , o* := — , Uj := Pj ■ Uj 

Pj Ci Ci 

Note that • x and • x are equal, but the order of the transformed values Cj may be 
different from the original order of the values Cj. A special case are the constraints 

i 

^ ■ Xj < bi, 

i=i 

with g > 0, by setting e* := g* and pj := q~fi These constraints can be used to model 
storage units including losses. 


Appendix B. Regression analyses 


In this section, we illustrate the regression analyses applied to interpolate the efficiencies 
in section [2] as well as the cost functions in section [31 

To calculate hourly energy loss (or corresponding energy reduction) factors for different 
storage capacities, the energy content of the storage unit is assumed to decrease exponen¬ 
tially. For simplicity, the thermal storage medium is assumed to be water with a specihc 
heat capacity of 4200 J/(kg-K) and the maximum possible temperature difference of the 
medium is set to 60 Kelvin. This implies a specihc energy capacity of 70 kWh/m^. Ad¬ 
ditionally, the energy is assumed to be homogeneously distributed within the storage. As 
a data basis for efficiencies of different storage units, the capacities (in liters) and the cor¬ 
responding daily losses (in kWh per day) were taken from specihcations of state-of-the-art 


products of a manufacturer of heat storage units j22|, pp. 24-25]. For the regression, these 
capacities were transformed from liter into kWh using the specihc energy capacity. The 
resulting data is listed in table IB.21 As a consequence of the observed linear correlation 
between the logarithmic values of the capacities and the corresponding daily energy losses, 
an exponential regression analysis based on the least squares method is used to model hourly 
energy reduction factors for diherent capacities. The result of the regression yields to the 
approximated hourly energy reduction factors gapprox; which are given by 


hour 

^approx 


0.2431954 24 


(B.l) 
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product 

1 

2 

3 

4 

5 

6 

7 

capacity [kWh] 

3.5 

5.6 

7.0 

8.4 

14.0 

21.0 

28.0 

daily losses [kWh/24h] 

0.54 

0.66 

0.79 

0.92 

1.4 

1.6 

1.8 

hourly reduction factor 

0.9930 

0.9948 

0.9950 

0.9952 

0.9956 

0.9967 

0.9972 

fitted factor 

0.9932 

0.9944 

0.9949 

0.9952 

0.9961 

0.9967 

0.9971 


Table B.2: Storage capacities and corresponding daily energy losses from a common manufac- 
tnrer of heat storage nnits and the corresponding honrly energy rednction factors compared 
with the resnlts from the exponential regression analysis, eqnation fIB.ip . 


where c is the capacity of the storage in kWh. This fnnction is nsed for a capacity-specihc 
consideration of losses in section |3l The data of the snpporting points from the prodnct 
specihcations and the resnlting htted honrly energy rednction factors are presented in ta¬ 
ble El 

The procednre for the estimation of the annnal costs for different charging and storage 
nnits nsed in sections 13.21 and 13.31 was as follows: The costs for the charge nnits stems from 
a linear regression analysis applied on the net valnes of the data from |23|, p. 26]. The resnlt 
is a linear fnnction with a slope of abont 7 [EUR/kW] and an offset of 130 EUR. For system 
optimization (i.e. minimization of (E])) the constant offset of 130 EUR was not considered, 
since it does not have an impact on the optimal conhgnration, which is dne to the linearity of 
the problem. Taking into acconnt a typical lifetime of 15 years the annnal costs for storage 
nnits, hence, were estimated at abont 0.47 EUR/kW. 

The cost-fnnction for storage nnits is the result of the averaged net values for heat storage 
units, neglecting installation and maintenance costs, which leads to 0.8 - 1.57 EUR/1 j^, 25 . 
For simplicity, the capacity costs are set to 1 EUR/1. In order to obtain annual costs per 
kWh, these costs were scaled with the storage unit specihc energy capacity and an assumed 
lifetime of 15 years, resulting in costs of about 0.95 EUR/kWh. 

Above considerations and data were used to define the specihc cost function (jl]) used in 
the examples from the more general function ([3]). 
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Appendix C. Python code 

#!/usr/biii/pythoii 

# coding: utf-8 

# (C) 2015 by Lars Siemer, NEXT ENERGY, EWE Research Center for 

# Energy Technology, 26129 Oldenburg, Germany. 

# Contact information: email: lars.siemer@next-energy.de 

# Optimization algorithm to solve the following constrained 

# linear program: 

# min c*x, 

# subject to: 0 <= x_i <= u_i 

# and a_i <= x_l+...+x_i <=b_i, for all i=l,...,n. 

# This code evolved in connection with the article 

# " Cost - Optimal Operation of Energy Storage Units: 

# Benefits of a Problem-Specific Approach" 
#================================================================ 

import numpy 

#================================================================ 

# definition of the algorithm 

def Algorithm(a,b,c,u): 

sigma=numpy.argsort(c) 
dim=numpy.size(c) 
x=numpy.zeros(dim) 
for k in range(0,dim): 
i = sigma [k] 

M_l=numpy.max(numpy.insert(a[:i] ,0,0)) 

M_2 = numpy.max(numpy.insert(a[i:] ,0,0)) 
m=numpy.min(b[i:]) 
if c [i] <0 : 

X [i] =numpy.min([u [i] ,m-M_l]) 
else : 

X[i]=numpy.min([numpy.max([0,M_2-M_l]),\ 
numpy.min([u[i],m-M_l])]) 
for 1 in range(i,dim): 
a[1] =a[1] -X [i] 
b[l]=b[l]-x[i] 

return x 
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